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Abstract. We present a review of cosmological solutions in non-linear massive 
gravity, focusing on the stability of perturbations. Although homogeneous and 
isotropic solutions have been found, these are now known to suffer from either Higuchi 
ghost or a new non-linear ghost instability. We discuss two approaches to alleviate 
this issue. By relaxing the symmetry of the background by e.g. breaking isotropy 
in the hidden sector, it is possible to accommodate a stable cosmological solution. 
Alternatively, extending the theory to allow for new dynamical degrees of freedom can 
also remove the conditions which lead to the instability. As examples for this case, 
we study the stability of self-accelerating solutions in the quasi-dilatonic extension 
and generic cosmological solutions in the varying mass extension. While the quasi- 
dilaton case turns out to be unstable, the varying mass case allows stable regimes of 
parameters. Viable self-accelerating solutions in the varying mass theory yet remain 
to be found. 



On the cosmology of massive gravity 2 

1. Introduction 

The search for a consistent theory of finite-range gravity is a longstanding problem and 
well motivated by both theoretical and observational considerations. On the theoretical 
side, whether there exists such a consistent extension of general relativity (GR) by 
a mass term is a basic question of classical field theory. After Fierz and Pauli's 
pioneering attempt in 1939 [1], this issue has been attracting a great deal of interest. 
On the observational side, continuing experimental probes of gravity have revealed new 
unexpected phenomena at large scales. One of the most profound discovery is the cosmic 
acceleration, which was found in 1998. The extremely tiny energy-scale associated with 
the cosmic acceleration hints that gravity might need to be modified in the infrared 
(IR). The massive gravity is one of the most interesting attempts in this direction. 

However, consistency of the massive gravity theory has been a challenging issue for 
several decades. Fierz and Pauli's model [1], which extends GR by a linear mass term, 
suffers from the van Dam-Veltman-Zakharov discontinuity [21 E]. Non-relativistic matter 
and relativistic matter couple to gravity with different relative strengths, no matter how 
small the graviton mass is. Although this problem can be alleviated by non-linearity, 
which was suggested by Vainshtein [1], the same non-linearity leads to a ghost instability. 
Indeed, at the non-linear level, the theory loses not only the momentum constraint but 
also the Hamiltonian constraint and, as a result, the non-linear theory includes up to 6 
degrees of freedom in the gravity sector. While 5 of them properly represent the degrees 
of freedom of a massive spin-2 field in a Poincare invariant background, the sixth one is 
the so called Boulware-Deser (BD) ghost [5]. 

Adopting the effective field theory approach in the decoupling limit, it was found 
that the BD ghost is related to the longitudinal mode of the Goldstone bosons associated 
with the broken general covariance [6]. Finding a theory free from the BD ghost was 
only recently achieved by de Rham, Gabadadze and Tolley (dRGT) [3 [8]. It was shown 
that the Hamiltonian constraint and the associated secondary constraint are restored in 
this theory and, as a result, the BD ghost mode is eliminated P [TOl [m[T2]. 

In order for the theory to be theoretically consistent and observationally viable, 
however, the absence of the BD ghost is not sufficient. At the very least, a stable 
cosmological solution is needed. 

The purpose of this article is to review the construction and the stability of 
cosmological solutions in the context of non-linear massive gravity. We start with 
describing the action of dRGT theory in Sec. [2l In Sec. [3] we construct homogeneous 
and isotropic cosmological solutions that exhibit self-acceleration. We then argue in 
Sec. m that all homogeneous and isotropic solutions in the dRGT theory are unstable 
and thus cannot describe the universe as we know it. In Sec. [5] we propose three 
alternative cosmological scenarios to avoid instabilities. One of them is based on the 
observation that breaking isotropy in the hidden sector (fiducial metric) still allows 
isotropic evolution of the visible sector (physical metric) and thus the standard thermal 
history of the universe. The other two proposals maintain isotropy in both visible and 
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hidden sectors but are based on extended theories of massive gravity with extra degrees 
of freedom, such as the quasi-dilaton theory [13] and the varying mass theory [T3] . 

2. Action 

We start by describing the action of dRGT massive gravity theory [8]. In order to have 
a manifestly diffeomorphism invariant description of the massive gravity, the action is 
built out of four Stiickelberg scalar fields, 0°(a;), a = 0, 1, 2, 3. These four scalars enter 
the gravity action through a "fiducial metric" defined as 

/m. = Umd^rdu<p' , (1) 

where fab{<P^) is a metric in the field space. The action for the gravity sector is a 
functional of the physical metric g^y and the fiducial metric /^;^. A necessary condition 
for the theory to be free from BD ghost is that the action in the decoupling limit vanish 
up to boundary terms when restricted to the longitudinal part of the Stiickelberg fields. 
With this requirement, the most general mass term without derivatives of g^v and f^v 
is constructed as 

Imass[g^lv, ftiu] = Mpirn^ / d'^x^/^ (£2 + "3 A + "4^), (2) 



with 



A = I m' - [/C^]) , C, = l {[ICf - 3 [JC] [JC^] + 2 [/C^]) 



£4 = ^ ([/C]^ - 6 [ICf [)C'] + 3 [JC'Y + 8 [/C] [/C^] - 6 [JC']) , (3) 

where a square bracket denotes trace operation and 

K = 5::-{y^^)\- (4) 

It was shown that the theory is free from BD ghost at the fully non-linear level even 
away from the decoupling hmit [HI [IDl [HI [12] ■ 

3. FLRW cosmological solution 

With the non-linear massive gravity theory free from BD ghost in hand, it is important 
to study its cosmological implications. In this section we thus construct Friedmann- 
Lemaitre-Robertson- Walker (FLRW) solutions. 

3.1. Open FLRW solution with Minkowski fiducial metric 

The original dRGT theory respects the Poincare symmetry in the field space and thus 
the fiducial metric is Minkowski, i.e. fab = Vat = diag(— 1, 1, 1, 1). In this subsection we 
thus review the FLRW solution with the Minkowski fiducial metric [15]. This is the first 
non-trivial FLRW solution in the context of dRGT massive gravity. 

In order to find FLRW cosmological solutions we should adopt an ansatz in which 
both g^i, and f^^, respect the FLRW symmetry. Since the tensor /^,^ is the puUback 
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of the Minkowski metric in the field space to the physical spacetime, such an ansatz 
would require a fiat, closed, or open FLRW coordinate system for the Minkowski line 
element. The Minkowski line element does not admit a closed chart, but it allows an 
open chart. Thus, while there is no closed FLRW solution, we may hope to find open 
FLRW solutions. A fiat FLRW solution, if it exists, is on the boundary between closed 
and open solutions but it was shown in [16] that such a solution does not exist. For 
these reasons, in the following we shall seek open FLRW solutions. 

Motivated by the coordinate transformation from Minkowski coordinates to Milne 
coordinates, we take the following ansatz for the four Stiickelberg scalars. 



/ = f{t)^l + \K\5,,x^x^ , <P' = v^/(t)x* , (5) 

with i^ < 0. This leads to the open FLRW form for the Minkowski fiducial metric, 

f,u = Va,d,rduct^' = - (/(t))' 5J<5° + \K\f{t)'n,^6l6i . (6) 

As for the physical metric, we adopt the general open {K < 0) FLRW ansatz as 
g^^dx^'dx" = - Nitfdt^ + n^dx'dx^ , 
n,,dx^dx^ = dx' + dy' + dz^ + ^}^^l+fy + '^f , (7) 

1 — K[x^ + y + -2 ) 

Here, x° = t, x^ = x, x^ = y, x^ = z; ^i,v = 0, ...,3 and i,j = 1,2,3. Without loss 
of generality, we assume that />0,/>0, a>0, A^>0. By substituting the above 
ansatz into the Einstein-Hilbert action plus the mass term ([2]), the gravity action up to 
a boundary term can be written as 

3d^a 



Ig = I d'^xVn 



3\K\Na -— + m (L2 + a^L^ + 0^14 



N 



(8) 



where 



L2 = 3a{a - ,/\K\f){2Na -fa- N^/\K\f) , 
L3 = (a - ^/\K\f)\ANa - 3fa - Ny^\f) , 
L4 = ia-./\K\ffiN-f). (9) 



In addition to the gravity action, we also consider a general matter content so that the 

total action is /tot = Ig + -^matter- 

Note that since the above ansatz fully respects the FLRW symmetry, the (Oi) 
components of the equations of motion for g^^ are trivially satisfied, thus the variation 
of the action with respect to N{t) and a{t) should correctly give all the non-zero 
components of the Einstein equation. On the other hand, because of the identity [18] 



' "^ ' "',o.r, (10) 



the number of independent equations of motion for the Stiickelberg scalars is one. 

Now let us take the variation of the action with respect to f{t), which contains all 
non-trivial information about the dynamics of the Stiickelberg scalars. It leads to 

(d - v^^) [(3 - 2X) + ^3(3 - X)(l -X) + a4(l - Xf] = 0, (11) 
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where X = ^^/\K\f /a. This equation has three solutions. The first one is a = ^/\K\N 
and corresponds to an empty open universe, i.e. the open FLRW chart of Minkowski 
spacetime. Thus this solution is not of our interest. The remaining two solutions are 



a _ 1 + 2a3 + ^4 ± a/1 + as + "i - ^4 ,-,„^ 

J - —?==^±, ^± = : • l-l-^j 

^y\K\ as + a^ 

Note that these two solutions are singular in the limit K -^ 0. This is consistent with 
the result in [16j, i.e. the non-existence of fiat FLRW cosmologies. On the other hand, 
with K < 0, hy taking the variation of the action with respect to N{t) and using f lT2|) . 
we obtain the following modified Friedmann equation. 

3H' + ^ = p^ + A^, H^^, (13) 



a 



Na' 



where pm is the energy density of the matter sector and 

,2 



rrin 



A- - ' 



(as + a^y 



1 + as ± ^/ 1 + as + a| - 0:4 



a, — 2a4 ± (1 + as) -\/ 1 + as + ao — a4 



(14) 



In this way the graviton mass manifests as the effective cosmological constant A±. When 
A-t > 0, the system exhibits self-acceleration. By taking the variation of the action with 
respect to a{t) we obtain a dynamical equation, which is consistent with the above 
modified Friedmann equation and the standard conservation equation for matter. 

3.2. Flat/closed/open FLRW solutions with general fiducial metric 

In appendix of Ref. p/Tj, the open FLRW solution was generalized to flat/closed/open 
FLRW solutions by considering a fiducial metric of the general FLRW type. (General 
fiducial metrics were first considered in [18] and the absence of BD ghost in this general 
setup was proved by [12] •) In this subsection we describe the general solutions. 

The most general fiducial metric consistent with fiat {K = 0), closed {K > 0) or 
open {K < 0) FLRW symmetries is 

where n and a are general functions of (p^, and Qij{ip'') is defined as in ([7]) with (x, y, z) 
replaced by (c/?^, yj^, yj^), and the curvature constant K is now either zero, positive or 
negative. Here, we have used the notation ip"" instead of 0" to make it clear that this 
form of the fiducial metric may be achieved from the original form ([T]) by non-trivial 
change of variables (as we have explicitly seen in the previous subsection). As for the 
physical metric, we adopt the ansatz (JTj) with an arbitrary sign for K. 

Similarly to the case in the previous subsection, the equation of motion for the 
Stiickelberg fields allows three branches of solutions. In the general setup at hand, the 
first branch is characterized by aH = aHf, where H = a/{Na) and Hf = a/{na) are 
Hubble expansion rates of the physical and fiducial metrics, respectively. Unfortunately, 
this branch would not allow nontrivial cosmologies since it does not evade Higuchi 
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bound [12] and thus linear perturbations around the corresponding solution [TBI EDI EI] 
include ghost in the cosmological history. Thus, we shall not consider this branch and 
restrict our attention to the other branches. 

The other two branches are characterized by a = Xj-a, where X± is the same as 
in (IT2I) . For these two branches, the metric equation of motion is exactly the same as 
flT3|) with flT^ . Surprisingly enough, the modified Friedmann equation (including the 
value of the effective cosmological constant induced by the graviton mass term) does 
not depend on the properties of the fiducial metric at all. When K± > 0, the system 
exhibits self-acceleration. 

4. New non-linear instability of FLRW solutions 

In the previous section we have constructed fiat, closed and open FLRW solutions in non- 
linear massive gravity with a general FLRW fiducial metric. The construction allows 
three branches of solutions. However, the first branch characterized by aH = aHf 
suffers from the Higuchi ghost at the level of linear perturbations and thus does not 
allow a non-trivial cosmological history. In this section we thus consider the other two 
branches of solutions characterized by a = X±a. These solutions evade the Higuchi 
ghost, but unfortunately we shall see that a new type of ghost instability shows up at 
non- linear level [22]. Based on this result, we shall argue that all homogeneous and 
isotropic FLRW solutions in the dRGT theory are unstable. We shall then propose 
alternative cosmological scenarios in the next section. 

4.I. Lin ear perturbation 

In this subsection, following p7], we shall investigate linear perturbations around a 
general flat/closed/open FLRW solution (characterized by a = X±a) with a general 
FLRW fiducial metric and an arbitrary matter content. We shall see that time kinetic 
terms for 3 among 5 gravity degrees of freedom always vanish at the level of the quadratic 
action, signaling for necessity of non-linear analysis. 

We first define perturbations of four Stiickelberg scalars through the exponential 
mapping, truncating at the second order, as 

y,-^ = x" + tt'^ + ^vr^'SfcTT" + Oin^) . (16) 

We then perturb the physical metric as 

goo = -N\t)[l + 2<f)], goi = N{t)a{t)(3i, g,, = a^ (t) [n^j + hi,] . (17) 

We suppose that tt"", (p, f3i, hi j = 0(e). The following gauge invariant variables can be 
constructed out of Stiickelberg and metric perturbations. 

0-^0-i-9,(iV7r°), /3- ^ /3, + ^ Avr° - ^TT,, 
JM a JM 

h^ = hij - DiTTj - Dj-Ki - 2NH7T^nij, (18) 

where Di is the spatial covariant derivative compatible with Qij. 
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In Sec. Elwe have seen that the mass term acts as an effective cosmological constant 
at the background level. Hence, we define 

where K± is specified in ([H]), and expand Imass instead of Imass- This greatly simplifies 
the perturbative expansion. As shown in [TTj, upon using the background equation of 
motion for the Stiickelberg fields but without using the background equation of motion 

for the physical metric, the quadratic part of Imass is simplified as 

r2 



i±s = ^ / d'xNa^VnMl^ [{h-f - f^Kl^ , (20) 



where 



X± is given by ([I2]), h"" = Vt'^hjj, h^ = Q^^W^hli, and Vt^^ is the inverse of n^j. This is 
manifestly gauge- invariant. 

What is important here is that the gauge- invariance of ImLss was shown without 
using the background equation of motion for the physical metric. This means that 
Imass is gauge-invariant for any matter content (as far as the matter action does not 
depend on the Stiickelberg fields so that the Stiickelberg equation of motion is derived 
solely from the graviton mass term) and that the remaining part /q^ = l[J^ — Imass 
of the total (gravity plus matter) quadratic action Il^{ is also gauge-invariant by itself. 
Hence, the remaining part Iq^ never depends on the Stiickelberg perturbations in any 
gauge choices. Another important point is that ImLss shown in (!20l) does not depend on 
(jf and /3f , and hence does not include time derivatives of Stiickelberg perturbations. 
Therefore, for any matter content, the dependence of the total quadratic action on the 
Stiickelberg perturbations is completely given by ( 120|) and time derivatives of Stiickelberg 
perturbations do not enter the quadratic action at all. This completes the proof of the 
statement that time kinetic terms for 3 among 5 gravity degrees of freedom always 
vanish at the level of the quadratic action. This proof holds for any matter content [IT] . 

The absence of quadratic kinetic terms for 3 gravity degrees of freedom shown in 
this subsection implies that the self-accelerating FLRW solutions evade the Higuchi 
bound [19] and thus are free from ghost at the linearized level even when the expansion 
rate is significantly higher than the graviton mass. At the same time, this signals for 
the necessity of a non-linear analysis. 

4-2. N on-linear perturbation 

In order to understand the physical content of the FLRW background in dRGT massive 
gravity, we need to investigate the reason why the kinetic term of one of the scalar 
modes and two of the vector ones have a vanishing kinetic term. We have assumed that 
this feature does not hold in general for this theory, but rather it is a consequence of 
the symmetries of the FLRW background. 
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Therefore it is convenient to study another background which has less symmetries 
than FLRW, but does lead FLRW in some limit. Probably, the simplest implementation 
of such a background is the axisymmetric Bianchi Type-I class of metrics, which can be 
written as 

ds^ = -N^{t)de + a2(t)[e^'^Wrfa;2 + e'^^'^'H.jdy'dy^] , (22) 

where here i,j G {2,3}, and y"^ = y,y^ = z. It is evident that this manifold is not 
isotropic, however as a approaches 0, the spacetime reduces to flat FLRW. Since we 
want to study the reason why there are missing kinetic terms for the perturbed FLRW 
fields, we will not consider the above metric as physical. Rather, for the moment being, 
we will merely use it as a tool to study non-linear perturbations on FLRW. In other 
words, linear perturbations on a anisotropic manifold will give information equivalent 
to non-linear perturbation theory on FLRW. The goal of this section is to show how 
the missing kinetic terms will depend on a, namely we will study their signs and the 
dispersion relations of the associated perturbation fields. 

According to the properties of the perturbation fields under a rotation in the y—z 
plane, we can decompose such fields into scalar (a.k.a. even) and vector (a.k.a. odd) 
modes. In particular we can write the metric for the even modes as 

rfsLn = - ^^(1 + 2$)(^t^ + 2aNdt[e^''d^xdx + e-^diBdy'] 
+ a^e^'^il + ilj)dx'^ + 2a^e''d^di/3dxdy' 

+ a^e-^''[6i,{l + r) + didjE]dyW , (23) 

whereas the metric for the odd modes reads 

dsl^^ = - N^dt^ + 2ae-"Nvidtdy' + 2a^e''d^Xidxdy' 

+ a'^e^^'dx^ + a^e-^^'iSij + d^ihj))dy'dy^ . (24) 

where the transverse condition holds, i.e. d'^Vi = = (9*Aj = = d^hi. At the same time, 
we also need to specify the form of the fiducial metric f^y. Since, we want to study the 
FLRW limit, we need the fiducial metric to possess already the FLRW symmetries from 
the beginning. The anisotropic physical metric corresponds to the simplest deviation 
from an overall FLRW symmetry for the whole system. Therefore we assume 

so that we need to give the perturbations of the Stiickelberg fields according to even/odd 
mode decomposition. In fact, the even- type perturbations of Stiickelberg fields read 

00 = t + 7r° , (j)^ = x + d^y , (j)' = y' + d'n . (26) 

For the odd modes sector, we consider instead 

<P' = t,<p' = x,(P' = f + n\ (27) 

where Avr* = 0. 
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It is possible to define gauge-invariant fields for even perturbations as follows 
$ =^-—d, ( ^ \ B = B + — T - ^^ " E 



r e-^'' ae"^" 






* ^ ^ 2a{H - E) iV ^* 

^^--'-JWJ^y ^.-^^-e-(/3-^i.), (28) 

where we have defined S = a/N. 

Then we can proceed to integrate out all the present auxiliary fields. In general, 
we can integrate out three modes, that is $, B and x- However, in the dRGT theory, 
it is possible to show that also the field f^^ can be integrated out. Therefore there are 
only three independent fields describing the even-mode perturbations, so we need to 
study the kinetic matrix of the three remaining fields, ip, P-^, and Et^. As o" — ?■ 0, the 
eigenvalues of the 3x3 kinetic matrix reduce to 

gp4 1 — r^ 2p]^ 

where we have introduced r = an/{aN), Mq^ = m^(l— r)X^[(l+2a3+Q;4) — («3+a4)X], 
X = a/a, pl = ki/iae^'^) ~ ki/a, pt = kT/{ae~'^) ^ kx/a, k'^ = ky + k^, and 
p^ = pi + Pt- The first and most important consideration is that k,2 and n^ have 
opposite sign. This property implies that a ghost will always be present in the even- 
modes sector, as the manifold approaches the FLRW limit. Furthermore, both k,2 and k,^ 
vanish in the exact FLRW case. One could wonder whether these modes, if their mass 
is finite - but non-zero, can be integrated out in this same FLRW limit. If the masses 
are heavy, the corresponding modes can be integrated out and the ghost is harmless 
in general. Otherwise, the ghost will be physical and the theory - at least on FLRW 
backgrounds - will not be consistent. 
In fact, we find 

^1 - P' + M^w , 

2-1 



24a 



(10p2 + ply - SpIpI - (2p2 + 3p2^) 



l-r2 



^3 - - ^2 + -TTT- (2pi + ^Pt) , (30) 

i2a 

with a;|w| < in general. Since there is not any mass-gap (i.e. there always exists 

some value of the momenta for which the frequency vanishes) for w^g, we conclude that 

the ghost is physical and cannot be integrated out from the Lagrangian. Therefore 

the FLRW background is not viable in the dRGT massive gravity theory. Notice that 

the first mode corresponds to the massive gravitational wave. Even though the theory 

succeeds in removing the BD ghost and giving the tensor mode a mass, it does not 

accept a stable FLRW solution. 
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We conclude this section by studying the odd modes. For these modes, a procedure 
similar to the one followed for the even modes leads to two independent fields. Therefore 
we confirm the expected presence of 5 dynamical degrees of freedom for this theory 
(3 even modes + 2 odd ones). The kinetic terms, and the frequencies for these two 
independent odd modes read 



Q PlPt _ Q Pr ^GW 2 _ 2 I /i^2 2 _ 2 

2p2 ' ^^2- 4^i_^2^ 



«i = o 9 ' '^s = -rr, ^ c^ , Wi = p + Mgw , ^2 = CoddP , (31) 



where c^^^ = (1 — r^)/(2cr). Therefore we find a massive tensor mode, and an healthy 
massless propagating mode (at speed Codd); provided that (1 — r)a > 0. 

5. Towards healthy massive cosmologies 

The appearance of the non-linear ghost shown in subsection 14.21 originates from the 
fact that quadratic kinetic terms exactly vanish: the kinetic terms show up at the 
cubic order and can become negative. The disappearance of kinetic terms at the 
quadratic order was shown upon using the background equation of motion for the 
Stiickelberg fields but without using other background equations. One can actually 
show that the off-shell quadratic kinetic terms have coefficients proportional to J<^ = 
(3 — 2X)+a3(3 — X)(l — X)+a4(l — X)^, where X = a /a, and that the self-accelerating 
FLRW solution is characterized by J^ = (or X = X± with X± shown in (TT2|) ). 

For this reason, in order to find a stable cosmological background, one needs to 
detune the proportionality between the quadratic kinetic terms and the Stiickelberg 
equation of motion characterizing the self-accelerating background. One way to achieve 
this would be to relax the FLRW symmetry by deformation of the background. 
This possibility will be considered in subsection 15.11 We shall find that relatively 
large deformation by anisotropy in the hidden sector (fiducial metric) may render 
the background solution stable. Another possibility would be to maintain the FLRW 
symmetry but to change the Stiickelberg equation of motion by adding extra dynamical 
degrees of freedom to the theory. We shall thus consider the quasi-dilaton extension in 
subsection 15. 21 and the varying mass extension in subsection 15.31 While self-accelerating 
FLRW solutions in the quasi-dilaton theory turn out to be unstable, the varying mass 
case allows some stable regimes of parameters. 

5.1. Anisotropic FLRW solution 

As argued above, the appearance of the non-linear ghost shown in subsection 14.21 is a 
consequence of the FLRW symmetry and the structure of the theory; in order to obtain 
a stable solution within the same theory, the FLRW symmetry needs to be relaxed. 

An inhomogeneous background solution was obtained in Ref. [16], where the 
observable universe is approximately FLRW for a horizon size smaller than the Compton 
length of graviton. Similar solutions with inhomogeneities in the Stiickelberg sector, 
meaning that the physical metric and the fiducial metric do not have common isometrics 
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acting transitively, were found in [23] . Note that those inhomogeneous solutions cannot 
be isotropic everywhere since isotropy at every point implies homogeneity. Note also that 
cosmological perturbations can in principle probe inhomogeneities in the Stiickelberg 
sector. For example, generic spherically-symmetric solutions are isotropic only when 
they are observed from the center of the universe. 

The goal of this subsection is, following |2l], to introduce an alternative option, 
where the assumption of isotropy is dropped but homogeneity, i.e. the cosmological 
principle, is kept. In a region with relatively large anisotropy, we find an attractor 
solution. On the attractor, the physical metric is still isotropic, and the background 
geometry is of FLRW type. Hence, the thermal history of the standard cosmology can 
be accommodated in this class of solutions. However, the Stiickelberg field configuration 
is anisotropic, which may lead to effects at the level of the perturbations, suppressed by 
smallness of the graviton mass. 

5.1.1. Fixed point solutions In this subsection, we review anisotropic FLRW solutions 
in the dRGT theory with the de Sitter fiducial metric [23]. The fiducial metric is 
obtained by taking (ITSll and setting fC = and Hf = a/ [no) = constant. For the 
physical metric we consider axisymmetric anisotropic extension of a fiat FLRW, i.e. the 
axisymmetric Bianchi type-I metric, given by 

g'^^^dx^'dx" = -dt^ + a\t) [e^'Wrfx^ + e-^''^''>6ijdfdy^] , (32) 

where indices i,j = 2, 3 correspond to the coordinates on the y—z plane. 
The Stiickelberg equation gives 

/f (if + 2S-i7/e"2"X) +2jf {H -^-Hfe"" X) = 0, (33) 

where 

jf^ = (3 + 3^3 + "4) - 2 (1 + 2a3 + 04) e"X + (03 + 04) e^'X^ , (34) 

J^^) = (3 + 3a3 + ^4) - (1 + 2^3 + a4)(e-'" + e")X + (03 + ai)e-''X^ , 
H = d/a, S = 0" and X = a/a. The metric equations of motion are given by 

3 {H^ - S^) - A = m^ [-(6 + 4^3 + 04) + (3 + 3 03 + a^) (2 e" + e-^'')X 

-(1 + 2 03 + 04) (e'" + 2 e-") X^ + (03 + 04) X^] , 

2 

t + SHE = -^{e-'^" - e'')X [(3 + 3 03 + 04) - (1 + 2 03 + a^{(f + r)X 



where 




+ (a3 + "4) re'X^j , (35) 



(36) 



We now look for solutions that are isotropic (S = 0, (T = ctq) and undergo a de Sitter 
{H = 0, H = Hq) expansion, which implies that the remaining parameters are also 
constant, i.e. X = Xq, r = tq. Excluding the fixed points which are isotropic ctq = and 
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those which exist only for special values of parameters (i.e. in a measure-zero subspace 
of the parameter space), we find a characteristic relation for the anisotropic fixed point, 
tq = e~^°"°, or, 

Xo = ^e"'\ (37) 

The remaining two equations allow us to determine Xq and ctq- One can show that this 
solution is stable against homogeneous linear perturbations if 



Ml = — . ^ ° J > , 3^ 

M4 + 9 1/2(^3 M2-M2) 

where [25] 

M^ = -^ [/7o'e'""(l + 2e=^'^«)(a3 + «4) -2HoHf{l + e^'^" + e^'^«) 



X 



;i + 2 as + a4) + ^J(2 + e3"")(3 + 3 as + "4)] 



+H]{3 + 3a3 + a4)] . (39) 

5.1.2. Linear perturbations The linear perturbation theory around axisymmetric 
Bianchi type-I backgrounds in dRGT theory was formulated in [25]. The formulation 
can be used to calculate the coefficients of the kinetic terms for the five gravity degrees 
of freedom. Similarly to the off-shell kinetic terms around isotropic FLRW solutions 
mentioned in the second paragraph of Sec. [5l it turns out that the kinetic term of one 
of the five degrees is proportional to Jv^ and other two are proportional to Jj^ . The 
remaining two correspond to the standard two polarizations of the tensor graviton and 
thus they always have finite and positive kinetic terms. 

By studying (!33|) . we see that on the fixed point (1371) . J^ 7^ while J^ = 0. 
Hence, the kinetic terms for two of the expected five gravity degrees of freedom vanish, 
signaling for necessity of non-linear analysis. 

5.1.3. Non-linear perturbations Due to the broken 5*0(3) symmetry, we can no 
longer use the standard scalar/vector/tensor decomposition for the perturbations. 
However, the axisymmetry of the background allows us to use 5*0(2) symmetry in 
the classification. The following analysis is based on Ref. [25] . 

Even modes. — The even mode perturbations are introduced according to the 
decomposition ( !23l) and once the non-dynamical degrees are integrated out, there are 
generically three dynamical degrees of freedom. However, once the background is fixed 
to be the anisotropic attractor solution, one of these modes has a vanishing kinetic term. 
On the other hand, we can still analyze the properties of the higher order kinetic terms 
by considering homogeneous deformations around the attractor solution, characterized 
by 

(T = ao + e(Ti + C(e2), S = eSi + ^(e^) = eai + ^(e') , (40) 



On the cosmology of massive gravity 



13 



where the background-physical- metric coefficients are given by gu = —N"^, gxx = o^e^'^, 
Qyy = 9zz = a^e'"^" ■ Furthermore, we have defined S as S = d, and expanded both a 
and S on the attractor solution. After diagonalization, the kinetic terms become 

-\ -1 



Kl 



^2 



f^Z 



Pt 



8M^ 



M4 + 9i72(3M2-M2) 
2 4 e^'^o M^pI [9 H^p'^ (M2 _ 3 m2) + M^pI (-2^^ + pi) 



(41) 



M^pKM^ - 3p2)2 _ 9H^(M^ - 3 M2) 6p^ + M^{-Ap'^ + pi) 
3 M2 e^'^" a^p^ [3 M2 (9 H^ + M^)ai + 2Ho{9M^-2 M^)T.i 



;i_e6'^o) M4_27i7oM3M2-M2) 



where p^ and px are the components of the physical momentum vector along the i. 
direction and on the y—z plane, respectively, while p^ = pi + p"^. Furthermore we have 
introduced the mass scale M^ as in Eq. f l39|) . Generically, the absence of ghosts imposes 
momentum dependent conditions. However, one can ensure stability at all scales by 
adopting the sufficient condition 

M^ < , M^ < - ^ ° 



27 Hi 



<0. 



(42) 



under which, both ni and k,2 can be made positive. For a parameter set which satisfies 
[2|) . the no-ghost condition for the third mode (with order e kinetic term) becomes 

-3\M^\ (9H^- \M'^\]ai + 2Ho (2\M'^\ -9|M2| 



O-Q 



Si 



<0 



(43) 



which depends linearly on the homogeneous deformations o"i and Si around the fixed 
point. Thus, regardless of the value of M^ and M^, there could always be a region where 
(Ti and Si conspire to render the third mode a ghost. On the other hand, if the initial 
conditions are such that the system is close to the attractor, it is possible to connect 
the evolution of Si algebraically to that of ai and obtain a regime where one can avoid 
the instability. By considering the equation of motion for o"i. 



Si + 3ifoSi + M2ai = 0, 



(44) 



we first note that the condition (!38|) for the stability of the fixed point against 
homogeneous perturbations, combined with the conditions ( H2l) . yields 

9Hi-\M'^\>0. (45) 

To satisfy the condition fH3l) . we suppose that the system is in the attractor regime, so 
that Si oc (Ti and that ai does not change sign during the course of evolution. This 
scenario can be attained if the friction term in f l44D dominates over the mass term, i.e. 

(46) 

Then, solving Eq. f H4|) and evaluating the solution at late times, we find the relation 

3 



9Hl>AMl 



Si 



H, + \l\Hl-Ml\a, 



(47) 



On the cosmology of massive gravity 14 

Thus, in this regime, the condition f l43p can in principle be satisfied by choosing the 
appropriate sign for ai. 

Odd m,odes. — The even mode perturbations are introduced according to the 
decomposition (p^ and once the non-dynamical degree is integrated out (another mode 
can be gauged away), there are generically two dynamical degrees of freedom. On the 
anisotropic attractor solution, the kinetic term of one of these modes vanishes, and as 
we did for the even modes, we consider homogeneous deviations from the fixed point to 
determine conditions for non-linear stability. 

After diagonalization, the kinetic terms become 



2p^ 
3 M2 e^'^o al p% [3 M2 (9 Hi + M'^)ai + 2Ho{9M'^-2 M^)^ 



4(l_e6'^o^ 



M^ - 27 Hi {3 M^ - M^) 



The first kinetic term is always positive, whereas the second mode acquires a kinetic 
term proportional to the deviation from the fixed point. In fact, up to a numerical 
factor, K2 above is the same as the kinetic term of the third mode in Eq. (H2|) . so if the 
conditions discussed in the even sector are satisfied, the odd sector will also be stable. 

5.2. Extended theory I: quasi-dilaton 

The simplest extension of dRGT theory is to introduce a scalar field a associated with a 
dilatation-like global symmetry (a has different meaning than the previous subsection) 

a^a- aMpi , 0" ^ e° 0" . (49) 

The action compatible with this symmetry is given in Einstein frame as |l3j I3. 

S=^ d^x v^ [R-2A + 2ml {C2{K:) + «3 C^iJC) + a^ C^iJC)) 

~ T?2" ^M^ d^ai- ^matter , (50) 

where £2, -^3 and £4 are given in Eq. (j3]), but the building block tensor (^ is replaced 
with 

/c:; -> £:; = Si: - e'^/^- [v^a^X ■ (51) 

5.2.1. Self- accelerating solutions We adopt the Minkowski fiducial metric and the flat 
FLRW ansatz for the physical metric as 

/^^ = -n'^(t)6l6'l + 6ij6l6i , g^^^dx" dx" = -dt^ + a\t) 6ijdx' dx^ . (52) 

The equations of motion for the Stiickelberg fields yield 

(l-X)X [3 + 3 ag + "4 - (3 as + 2 ^4) X + a4 X 2] = , (53) 

a^ 

I There is an additional term allowed by the symmetry, J d'^x -y/— / e^ '^/^v ^ which does not change the 
conclusions of the present discussion (See [1S1I27] for details). 
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where X = e'^/^^^/a, leading to four attractors: X = 0, X = 1 and X = X± with 



X, 



3 tta + 2 ^4 ± a/9 a| — 12 a^ 
2 0:4 



(54) 



Among them, X = and X = 1 leads to either strong coupling or instability [T3] . 
We thus consider X = X± as backgrounds. Along these branches of solutions, the 
(modified) Friedmann equation becomes 



3-|)//^ = A + A±, 

where the graviton mass manifests as the effective cosmological constant 

,2 



A 



m„ 



± 



2a| 



[9 {3aj-6alai + 2 aj) ± a^ (9 al - 12 a^f^^] . 



(55) 



(56) 



From (l55!l . we immediately see that a sensible cosmology requires u < 6. Finally, the 
equation of motion for the quasi-dilaton field gives 



na 



1 + 



m2X2[a3(X-l)-2]- 



(57) 



5.2.2. Perturbations We now introduce perturbations as [26] 
= Mpi [log(a X) + 5a] , (5^oo = -2 $ , (5c/oi 



0" 



^9i, 



a{BI + diB) 



2 ^ 5,, + ( 9.9, - ^ 9,9') i? + \{d,Ej + 9,En + /.J^ 



, (58) 



0, where Bj and Ej are transverse and hj^ 



3 / 2 

while we fix the unitary gauge 
is transverse and traceless. With respect to the dRGT theory, we have an additional 
scalar field, so in total, we expect 2 tensor, 2 vector and 2 scalar degrees, once the 
non-dynamical modes are integrated out. 

Tensor perturbations. — The quadratic action for the tensor modes reduces to 



S-T 



Ml, 



(Tk a^ dt 



\h. 



TT\2 
ij I 



^ + M, 



GW 



\h 



TT\2 
ij I 



where 



m2(r-l)X3 



+ 



(59) 



(60) 



X-1 \r-l X-1, 

Generically, Mqw ~ ^{H) so even if the tensor modes are tachyonic, the time scale of 
their instability is of the order of the age of the universe p6] . 

Vector perturbations. — For the vector modes, the quadratic action is 



S^ 



Ml, 
16 



(fk a'^ dt k 



IK 



T|2 



1 + 



fc2(r2-l) 
2a2 _ff2 ^ 



-M^^lEfP 



(61) 



We see from (!6T!) that if (r^ — l)/uj < 0, there is a critical momentum above which the 
vector modes have ghost instability. Therefore, the UV cutoff scale of the effective theory 
Kuv should be lower than this critical (physical) momentum to ensure the stability of 
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the system: A^^ (1 — r^)/{H'^uj) < 2. In addition, the frequency in the canonical 
normahzation yields a further condition on avoiding tachyonic instability, which arises 
if Mqy^^ > and (r^ — l)/u) > 0. The growth rate of this instability can be made lower 
than or at most of the cosmological scale for all physical momenta below the UV cut-off 
Auv, provided that [26] 

2H^u 

Scalar perturbations. — After integrating out 6gof^ as well as the would-be Boulware- 
Deser degree of freedom, the scalar sector contains two coupled modes. The kinetic part 
of the quadratic action is formally written as 



^Iv^zF^- (62) 



Ss 3 / ^^-^ a^ dt 



KuYl' + K22Y2' + Ki2 YiY2 +^2^1 



, V * 



(63) 



where Yi and Y2 are particular linear combinations of \E' and 6a. For our purposes, it is 
enough to study the determinant of the kinetic matrix, given by 

3 fc6 a;2 a^ H^ 



det K = KnK22 - Kl^ = t "'^ " ^ T , (64) 



ua^H^-4^, 

(6-w) 



r-1 



The absence of ghost degrees in the scalar sector requires det i^ > as a necessary 
condition. We first note that the determinant is always negative if w < 0. Along with 
the condition obtained from ( 155|) . we thus obtain < w < 6 as a necessary condition. 

Furthermore, demanding that det K > {) for all physical momenta below the UV 
cutoff of the theory, we obtain 



^< 2 <2' ('') 

where we have used the condition < a; < 6 to obtain the last inequality. Unfortunately, 
(1651) is not acceptable since it would imply that the UV cutoff scale would be lower 
than the cosmological scale and that the theory would not be applicable to cosmology. 
Therefore, we conclude that for physical wavelengths shorter than cosmological scales, 
det K <Q and one of the two degrees of freedom is a ghost [26] . 

It can also be checked (see [26] for details) that energies of the ghost mode are not 
parametrically higher than H ~ mg. This signals the presence of ghost instabilities in 
the regime of validity of the effective field theory. 

5.3. Extended theory II: varying mass 

A further way of extending the dRGT theory is to allow the parameters of the theory 
to vary with a scalar field a. The action in this case is [1] 



S= j d^x v^ Jm^i [^ - A + mlia) [A + «3(^) A + a^{(j) £4] 



-]^d^ad''a-V{a) + Crn\ , (66) 



where £2, £3 and £4 are given by Eq. 
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5.3.1. Background As in the previous subsection, we adopt the Minkowski fiducial 
metric and the fiat FLRW ansatz for the physical metric fl52l) . The equations of motion 
for the Stiickelberg fields yield 

<{X - 1) 



[3 — 3 (X — 1) 0:3 + (X — 1)^ 04] = constant, 



(67) 



X3 

where X = 1/a and r = an. Due to the assumptions of fiat space and Minkowski 
fiducial metric, if rrig and 03^4 are time- independent, the solution X = constant does 
not allow any non-trivial cosmologies (see the second paragraph of subsection I3.ip . 
By defining 

p„ = M^yml (X - 1) [6 + 4 03 + 04 - X (3 + 5 ^3 + 2 04) + X^ (03 + a^ 



Pn 



Q 



Ml^m] [6 + 4 03 + 04 - (2 + r) X (3 + 3 «3 + 04) 



M|i ml b 



+ 



(T 



2 m' 



rUn 



+ (1 + 2 r) X^ (1 + 2 03 + 04) - r X^ (03 + a^ 
(X - l)2|a3 (4 - X - 3rX) + a'^{X - 1) (rX - 1) 

3 _ (X - 1)«3 + '^^^ [3 - 3 (X - 1) «3 + (X - If a,] I 



P.-y + V. 



Pa 



2 



(6J 



we can write the set of background equations of motion in the following form 

I , . .V 1 



3H' = A 



Ml, 



{pa + Pr. 



H 



2 Mil 



[{Pa+Pa) + {Pm+Pm)] 



Pm + SH{pm+ Pn 



-Q, p^ + 3H{p„+p^) = Q. 



(69) 



where prime denotes differentiation with respect to a. 

Although dynamical analysis for these equations have been studied in the literature 
[2S1 [2S], there is not yet a simple self-accelerating solution in the varying parameter 
massive gravity. In the following, we do not assume any specific evolution and keep the 
functions mg{a), 03(0") and 04(0") generic. 



5.3.2. Perturbations We now introduce perturbations, following [SB]. The metric is 
decomposed as in Eq. (l58ll and we adopt the unitary gauge as (50° = 0, while the scalar 
field is perturbed as 

a= {a) + Mpi 5a . (70) 

Tensor perturbations. — The tensor action reduces to 



s^ 



Ml, 



d'^k a'^ dt 



\h. 



TT\2 



+ M, 



GW 



2 \ \vTT\2 



where 



Mlw 



(r-l)X^ 



^l{x-i] 



Ml,, 



V 



1 2X 

+ 



(71) 



r-1 X- 1 



Pm ~r Pn 

Ml, 



.(72) 



The stability of long wavelength tensor modes is ensured by M^^ > 0. 
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Vector perturbations. — For the vector modes, the action is 



S^ 



V 



Mh 


f d^ka^dtk^ 




IBP 


-M'aw\El\' 


16 . 


l(- 


fc2(r2_l)M2j\ 

2a2 (pm+Pm) J 



(73) 



By requiring that the kinetic term is positive for all physical momenta below the cut-off 
scale of the theory Ajjv, we obtain the stability condition for the vector modes as 



A2 



H^R 



<2, R 



Pm ~r Pm 



(74) 



Under the above condition, we can further analyze the stability of the vector sector, 
by introducing a time reparametrization which renders the modes Ef canonical, then 
requiring that their frequency is an increasing function. This procedure yields sufficient 
(but not necessary) condition for stability 



1 + 



d , fRM: 
m 



GW 



8NH dt \r^ -I 



A?,v(l-r2) 



3 

^ 2 ^ ANH 



1 rfln(M2^) 



.(75) 



H^R ~ 2 ' ANH dt 

Scalar perturbations. — As in the quasi-dilaton theory, we integrate out the 
nondynamical degrees and are left with two coupled modes in the scalar sector. The 
kinetic part of the action is formally 



Ss 3 I -^ a^ dt 



KulYA^ + K22\Y2\^ + Ki2 



VlFa' + V-aFi' 



(76) 



where Yi and Y2 are linear combinations of \1' and 6a. For our purposes, it is enough to 
study det K = K11K22 — -ft'12, whose explicit form is 



detK 



3 M2j a^ k^ {pm + Pmf {pa + Pa 



6M^,H^] 



6M^^H^ 



(77) 



^ 1)2 [AM^.H^ g - (p^+p^)(p^+p^ 

By requiring that the determinant is positive, we see that in order to avoid a ghost 
degree of freedom, the momenta in the range < k/a < Auv should all satisfy 



Pa+Pa 

AMl^m 



1^ 



pm \ Pn 



(7J 



Explicit diagonalization of the system shows that this condition is actually a sufficient 
condition to avoid ghost instabilities in the scalar sector [26] . 

For a background solution which can effectively describe the late time acceleration, 
we can assume a de Sitter like expansion, i.e. \H\ -C H^. With these considerations, the 
stability requirement for the scalar sector becomes even simpler, 

4 e 



R + 



>0, 



(79) 



R-Q H^a^ 

where R is defined in f in|) . If indeed all the physical momenta below the cutoff scale 
Auv satisfy (179!) and if we suppose Auy/H > 3/2 so that the theory is applicable to 
cosmological scales, then the no-ghost condition for scalar perturbations in the regime 
\H\ <^ H^ becomes simply 



R>Q. 



(80) 
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6. Summary and Discussion 

The extension of GR by a mass term has been studied for several decades. Nonetheless, 
a self-consistent non-linear massive gravity theory with five propagating degrees of 
freedom, dubbed the dRGT theory, has been proposed only recently. 

In the present article, we reviewed several cosmological solutions in the context 
of the dRGT theory [71 [8]. We have firstly described open FLRW solutions with a 
Minkowski fiducial metric. By considering a general FLRW-form fiducial metric, the 
branch of open FLRW solutions was generalized to FLRW solutions with general spatial 
curvature. However, all these FLRW-type cosmological solutions lead the kinetic terms 
of three among five gravity degrees of freedom to vanish at the level of the quadratic 
action. This phenomenon is a consequence of the symmetry of the FLRW background. 
On analyzing the behavior of the non-linear perturbations, it was then shown that there 
is always at least one ghost (among the five degrees of freedom) in the gravity sector. 

We have then discussed two approaches towards healthy cosmologies in massive 
gravity. One proposal is to introduce relatively large anisotropy in the configuration of 
Stiickelberg fields, which form the hidden sector of the theory. We considered the fixed 
point solution named as "anisotropic FLRW" , a solution with the FLRW symmetry in 
the visible sector (physical metric) but with anisotropy in the hidden sector. Performing 
a non-linear analysis around the anisotropic fixed point yields that anisotropic FLRW 
solutions can be ghost-free for a range of parameters and initial conditions. The second 
proposal discussed here consists of introducing an extra degree of freedom coupled to 
the hidden sector. As examples for this possibility, we have considered the quasi-dilaton 
theory and the varying mass model. For the quasi-dilaton theory, the self-accelerating 
theory turns out to be unstable. On the other hand, in the varying mass case, there 
is a regime of parameters in which a stable cosmological evolution is possible, although 
viable self-accelerating solutions yet remain to be found. 

Besides the stability investigation, the study of observational signals from graviton 
mass, although not included in this review, is also important. For example, in [30], it 
was found that graviton mass may leave a prominent feature with a sharp peak in the 
stochastic gravitational wave spectrum. The position and height of the peak may tell us 
information about the graviton mass today and the duration of the infiationary period. 

Last but not least, as a developing field, massive gravity still leaves many intriguing 
unsolved questions. One of the most interesting questions is the construction of possible 
UV completion of massive gravity. One of the possible directions to this end would be 
to seek a mechanism that realizes the specific structure of the graviton mass term as a 
consequence of a spontaneous symmetry breaking. Another important question is the 
fate of super-luminal mode [31] in the gravity sector. It is generically expected that 
in the massless limit, observable effects of the super-luminal mode should disappear 
and GR should be recovered, provided that the mode is excited by a fixed amount of 
matter source. Thus, it should be possible to obtain an observational upper bound on 
the graviton mass although it is probably not stronger than m^ < 0{Hq). 
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